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ABSTRACT 

We present measurements of the velocity dispersion profile (VDP) for galaxy groups in the final data 
release of the Sloan Digital Sky Survey (SDSS). For groups of given mass we estimate the redshift-space 
cross-correlation function (CCF) with respect to a reference galaxy sample, £^(r p ,7r), the projected 
CCF, w p (r p ), and the real-space CCF, £ C g{r)- The VDP is then extracted from the redshift distortion 
in (r p , n), by comparing (r p , it) with £ C g(0- We find that the velocity dispersion (VD) within 
virial radius (-R200) shows a roughly flat profile, with a slight increase at radii below ~ 0.3i?2oo for high 
mass systems. The average VD within the virial radius, a v , is a strongly increasing function of central 
galaxy mass. We apply the same methodology to Y-body simulations with the concordance A cold 
dark matter cosmology but different values of the density fluctuation parameter eg , and we compare 
the results to the SDSS results. We show that the a v — M* relation from the data provides stringent 
constraints on both er g and a ms , the dispersion in logM* of central galaxies at fixed halo mass. Our 
best-fitting model suggests as — 0.86 ± 0.03 and a ms = 0.16 ± 0.03. The slightly higher value of as 
compared to the WMAP7 result might be due to a smaller matter density parameter assumed in our 
simulations. Our VD measurements also provide a direct measure of the dark matter halo mass for 
central galaxies of different luminosities and masses, in good agreement with the results obtained by 
Mandelbaum et al. (2006) from stacking the gravitational lensing signals of the SDSS galaxies. 
Subject headings: dark matter - galaxies: halos - large-scale structure - method: statistical 



1. INTRODUCTION 

Satellite galaxies are an important tracer of the po- 
tential well within which they reside, as they can be 
observed to much larger radii than other tracers. The 
potential well is determined by both the dark matter 
particles and baryons. In the past one and a half decade, 
Y-body simulations provided excellent understanding of 
the mass profiles and shapes of dark matter haloes. For 
example, most haloes are well-fitted by the Navarro et al. 
(1997, hereafter NFW) profile, while their shapes are 
triaxial ellipsoids (Jing & Suto 2002). However, it is 
expected that baryonic condensation at the centers of 
dark matter haloes may modify the inner profiles signif- 
icantly, particularly on the galaxy scale (White & Rees 
1978; Blumenthal et al. 1986). Satellite galaxies provide 
a valuable probe of the dark matter halo profiles and the 
baryonic condensation processes. 

On galaxy scale, satellites in the Milky Way provides 
the best dynamical constraint on its mass profile out 
to about 200 kpc (Kochanck 1996; Wilkinson & Evans 
1999). Stellar-dynamical analysis of massive early-type 
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galaxies also provide the kinematics of the inner few 
kpc, which seems to show the density profile is roughly 
isothermal (Koopmans et al. 2006). This is much steeper 
than the inner slope expected from the NFW profile, 
which suggests that baryonic processes have modified 
the inner profiles, consistent with theoretical expecta- 
tions (e.g. Gnedin et al. 2004; Lin et al. 2006). For clus- 
ters of galaxies, using the CNOC1 survey Carlberg et al. 
(1997b) showed that for their clusters, the line of sight ve- 
locity rises from 0.1 virial radius, reaches a peak around 
0.3 virial radius and then shows a roughly flat profile with 
a very slight decline. The line of sight velocity dispersion 
profile is consistent with that derived by Lokas & Mamon 
(2003) from a detailed study of the nearby Coma cluster. 
The CNOC1 data have been used by van der Marel et al. 
(2000) to derive the mass distribution of clusters. 

There have also been many studies that measure the 
average linc-of-sight (los) velocity dispersion of the satel- 
lite galaxies within groups or clusters in order to infer the 
mass of their host dark matter halo. Early studies were 
mostly limited to rich galaxy clusters with a large num- 
ber of satellite galaxies (e.g. Carlberg et al. 1996, 1997a). 
For less massive systems, halo mass estimates are usu- 
ally obtained for central galaxies of similar luminosity or 
stellar mass by stacking the kinematics of their satellite 
galaxies (e.g. Erickson et al. 1987; Zaritsky et al. 1993; 
Zaritsky & White 1994; Zaritsky et al. 1997). Taking 
advantage of the large redshift surveys, in particular 
the Two-degree Field Galaxy Redshift Survey (2dFGRS; 
Colless et al. 2001) and the Sloan Digital Sky Survey 
(SDSS; York et al. 2000), recent studies have applied 
this technique to large samples of satellite galaxies and 
studied the dependence of the los velocity dispersion on 
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galaxy properties such as luminosity, stellar mass and op- 
tical color (e.g. McKay et al. 2002; Brainerd & Specian 
2003; Prada et al. 2003; van den Bosch et al. 2004; 
Becker et al. 2007; Conroy et al. 2007b; Norberg et al. 
2008; More et al. 2009b, 2011). These studies have well 
established that the los velocity dispersion increases with 
galaxy luminosity and mass, consistent with the theoret- 
ical expectation that more massive galaxies are hosted 
by more massive halos. In addition, More et al. (2009b) 
found that the luminosity-halo mass relation may differ 
with color, while the stellar mass-halo mass relation is 
less color-dependent. 

In this work we measure the velocity dispersion profile 
of satellite galaxies in groups by modelling the redshift 
distortions in the two-point cross-correlation function be- 
tween galaxies and groups of galaxies. It has been well 
established that, the two-point autocorrelation function 
(2PCF) of galaxies measured from redshift surveys is dis- 
torted along the line of sight due to the peculiar motions 
of galaxies. On small scales the 2PCF is stretched along 
the line of sight, called 'Finger-of-God' (FoG) effect, and 
on large scales it is squashed due to the global infall 
of galaxies towards high-density regions (Kaiser 1987). 
Thus the redshift distortion (RSD) in 2PCF contains 
useful information about the relative motions of galax- 
ies, and has been used to measure the so-called pairwisc 
velocity dispersion (PVD) of different classes of galaxies 
by many authors (Davis & Peebles 1983; Mo et al. 1993; 
Fisher et al. 1994; Zurek et al. 1994; Marzke et al. 1995; 
Somerville et al. 1997; Jing et al. 1998; Zchavi et al. 
2002; Hawkins et al. 2003; Jing & Borner 2004; Li et al. 
2006a). These measurements have been used to 
test/constrain both semi-analytic models of galaxy 
formation (e.g. Li et al. 2007) and halo occupa- 
tion distribution models of galaxy distribution (e.g. 
Jing & Borner 2004; Slosar et al. 2006; Tinker et al. 
2007; van den Bosch et al. 2007), as well as cosmologi- 
cal parameters (e.g. Jing et al. 1998; Yang et al. 2004; 
Wang et al. 2008; Cabre & Gaztanaga 2009). 

Here we use the same methodology as used in these pre- 
vious studies, except that we use the galaxy-group cross- 
correlation function to probe the velocity distributions of 
galaxies with respect to dark matter haloes. When com- 
pared to the aforementioned 'direct' measurements, our 
method has the advantage that it is not influenced by 
the so-called 'interlopers', i.e. mis-identified group mem- 
bers (see Yang et al. 2005a, for detailed discussion) . The 
cross-correlation functions between groups and galaxies 
have been determined by Yang et al. (2005b) based on 
an earlier data release of the Sloan Digital Sky Survey 
(SDSS; York et al. 2000). In this study we make use of a 
sample of about 16,000 groups which is constructed from 
the final data release (DR7; Abazajian et al. 2009) of the 
SDSS using the halo-based group finding algorithm of 
Yang et al. (2007). We estimate both the redshift-space 
cross-correlation function £,^(r p ,Tr) and the real-space 
cross-correlation function ^ C g{f), between a given sub- 
sample of groups and a reference sample of galaxies se- 
lected from the SDSS. The velocity dispersion profile for 
each subsamplc is then derived by modelling the redshift 
distortion in t;^ (r p ,Tr). As we will show, the velocity 
dispersion profiles within the virial radius can be reliably 
determined using the current data over a wide range in 



central galaxy luminosity and mass. More interestingly, 
comparisons with a set of high-resolution A-body simu- 
lations show that the average velocity dispersion within 
the virial radius as a function of central galaxy mass 
provides stringent constraints on both the density fluc- 
tuation parameter a$, and the correlation of the stellar 
mass of galaxies with the dark matter mass of their host 
halos. 

Throughout this paper we assume a cosmology model 
with the density parameter f2 m = 0.27 and the cos- 
mological constant fl\ = 0.73, and a Hubble constant 
H = lOO/ikms^Mpc" 1 with h = 0.7. 

2. DATA 
2.1. Group catalog 

The galaxy group catalog used in this paper is con- 
structed by Yang et al. (2007) from sample dr72 of the 
New York University Value-Added Galaxy Catalogue 
(NYU-VAGC). The NYU-VAGC is a catalogue of lo- 
cal galaxies (mostly below z ~ 0.3) selected from the 
SDSS data release 7 (DR7; Abazajian et al. 2009), pub- 
licly available at http://sdss.physics.nyu.edu/vagc/, and 
is described in detail in Blanton et al. (2005). The main 
virtue of the NYU-VAGC is that it provides a detailed 
account of the selection effects in the survey, thus suit- 
able for statistical studies of the galaxy distribution and 
large-scale structure in the local Universe. 

To select the groups of galaxies, a modified version 
of the halo-based group-finding algorithm developed in 
Yang et al. (2005a) is applied to a sample of ~ 6.4 x 10 5 
galaxies selected from the NYU-VAGC sample dr72 
with redshifts in the range 0.01 < z < 0.20 and with 
a redshift completeness above 70%. The reader is re- 
ferred to Yang et al. (2005a) and Yang et al. (2007) for 
detailed description of the group finder. The group cat- 
alogue contains about half a million systems, of which 
the majority have only a single member. In this work we 
use a subset of ^16,000 groups that have at least three 
member galaxies. 

We use the most massive galaxy member of each 
group as the group center which is called the "central" 
galaxy in what follows. The stellar mass of each cen- 
tral galaxy accompanies the NYU-VAGC release, which 
is estimated based on its redshift and the five-band mag- 
nitudes from SDSS photometric data, as described in 
detail in Blanton & Roweis (2007). This estimate cor- 
rects implicitly for dust and assumes a universal stel- 
lar initial mass function (IMF) of Chabricr (2003) form. 
As demonstrated in Appendix A of Li & White (2009), 
once all estimates are adapted to assume the same 
IMF, the Blanton & Roweis masses agree quite well 
with those obtained from the simple, single-color esti- 
mator of Bell et al. (2003) and also with those derived 
by Kauffmann et al. (2003) from a combination of SDSS 
photometry and spectroscopy. In this work we use the 
"total masses" instead of the "Petrosian masses" used 
by Li & White (2009), obtained by correcting the lat- 
ter using SDSS "model magnitudes" (see Appendix A of 
Guo et al. 2010, for details). 

From the group catalogue, we select 18 subsamples ac- 
cording to the stellar mass of the central galaxies, rang- 
ing from \og{M tt /Mo) =10.0 to 12.0. Each subsam- 
ple includes central galaxies in a stellar mass interval 
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TABLE 1 

SUBSAMPLES OF GROUPS SELECTED ACCORDING TO THE STELLAR MASS OF THEIR 
CENTRAL GALAXIES, INCLUDING THE STELLAR MASS RANGE, THE NUMBER OF GROUPS, 
THE MEAN STELLAR MASS AND THE MEAN LUMINOSITY OF EACH SUBSAMPLE, AS WELL 
AS THE ONE-DIMENSIONAL VELOCITY DISPERSION AND HALO MASS DETERMINED FOR 

EACH SUBSAMPLE. 



g(M,/M Q ) 


N group 


log«M*/M Q » 


log«L/L » 


a v /kms 1 


log(M ft /M ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


10.3, 10.6) 


913 


10.49 


10.33 


133±11 


12.17±0.13 


10.4, 10.7) 


1415 


10.59 


10.40 


139±8 


12.24±0.09 


10.5, 10.8) 


2190 


10.69 


10.47 


151±6 


12.36±0.06 


10.6, 10.9) 


3145 


10.78 


10.53 


166±5 


12.50±0.05 


10.7, 11.0) 


4245 


10.88 


10.61 


188±6 


12.69±0.05 


10.8, 11.1) 


5257 


10.97 


10.69 


216±7 


12.88±0.04 


10.9, 11.2) 


6214 


11.07 


10.77 


243±6 


13.04±0.04 


11.0, 11.3) 


6537 


11.16 


10.85 


291±8 


13.29±0.03 


11.1, 11.4) 


6014 


11.24 


10.93 


338±9 


13.49±0.03 


11.2, 11.5) 


4685 


11.33 


11.01 


385±10 


13.65±0.03 


11.3, 11.6) 


3115 


11.42 


11.09 


439±15 


13.82±0.04 


11.4, 11.7) 


1806 


11.51 


11.18 


522±22 


14.04±0.05 


11.5, 11.8) 


903 


11.60 


11.26 


591±32 


14.20±0.07 


11.6, 11.9) 


374 


11.69 


11.35 


654±38 


14.32±0.07 



of 0.3 dex, with successive subsamples overlapping by 
0.2 dex. As we will show below, the velocity dispersion 
profile cannot be reliably determined for subsamples be- 
low log(M*/M Q ) = 10.3. We will also ignore the most 
massive subsamplc with 11.9 < log(A/»/A/ ) < 12.0 
when measuring the velocity dispersion profile, because 
its redshift-space cross-correlation function is noisy due 
to the small sample size. The stellar mass range, the 
number of groups, and the mean stellar mass and lumi- 
nosity of the rest 14 subsamples are listed in the first 
four columns in Table 1 . The luminosity of each galaxy 
is computed from its r— band absolute magnitude and 
the absolute magnitude of the Sun is assumed to be 4.76 
mag following Blanton et al. (2003). 

2.2. Reference galaxy sample and random sample 

We have constructed a magnitude-limited galaxy sam- 
ple from the NYU-VAGC sample dr72, which will serve 
as our reference sample when measuring the group- 
galaxy cross-correlation functions. This consists of about 
half a million galaxies with r < 17.6, —24 < Mo.i r < —16 
and redshifts in the range 0.01 < z < 0.2. Here, r 
is the r-band Petrosian apparent magnitude, corrected 
for Galactic extinction, and Mo.i r is the r-band Pet- 
rosian absolute magnitude, corrected for evolution and 
X-corrected to its value at z = 0.1. These selection crite- 
ria, with the exception of the redshift range, are the same 
as in our previous papers where we studied the cluster- 
ing of galaxy luminosity and stellar mass (e.g. Li et al. 
2012). 

To obtain reliable estimates of the cross-correlation 
functions, which will form the basis of our study, the 
cross pair counts between the groups and the reference 
galaxies must be compared with the cross pair counts 
between the same set of groups and a 'random sample'. 
The random sample is unclustcrcd but fills the same re- 
gion of the sky and has the same position- and redshift - 
dependent selection effects as the reference sample. We 
have constructed our random sample from the observed 
reference sample itself, as described in detail in Li et al. 
(2006b). For each real galaxy we generate 10 sky posi- 
tions at random within the mask of sample dr72, and 



we assign to each of them the redshift of the real galaxy. 
The resulting random sample is valid for clustering anal- 
yses, provided that the survey area is large enough so 
that structures in the real sample are wiped out by ran- 
domizing in angle, and that the effective depth of the sur- 
vey does not vary from region to region. Both are true 
to good accuracy for our sample, which covers > 6000 
deg , is complete down to r = 17.6 and is little affected 
by foreground dust over the entire survey region. Exten- 
sive tests show that random samples constructed in this 
way produce indistinguishable results from those using 
the traditional method (Li et al. 2006b). 

2.3. Simulations and dark matter subhalo catalogs 

In this paper, we use three iV-body simulations and 
dark matter halo catalogs constructed from them to 
study the connections of galaxy velocity dispersions with 
halo masses and cosmological parameters. The simula- 
tions use 1024 3 particles to follow the dark-matter distri- 
bution in a cubic region with 300/i Mpc on a side, cor- 
responding to a particle mass of 1.87 x 10 9 /i _1 M Q . These 
are obtained using an upgraded version of the particle- 
particle-particle-mesh (P 3 M) code of Jing & Suto (1998, 
2002), which has incorporated the multiple- level P 3 M 
gravity solver for high-density regions (Jing & Suto 2000; 
Jing et al. 2007). The cosmology model assumed is a spa- 
tially flat ACDM model with the same parameters Vl m = 
0.268, fl A = 0.732, n b = 0.045, H = 71kms- 1 Mp C - 1 
and n = 1, except that a different density fluctuation pa- 
rameter is adopted, which is as = 0.75, 0.85 and 0.95, for 
the three simulations respectively. The simulations are 
described in detail in Jing et al. (2007). The cosmology 
for the simulations is nearly identical to the cosmology 
we assume in this study (mentioned in the introduction). 

Dark matter halos are identified for each of the simula- 
tion outputs using the standard Friends-of- Friends (FoF) 
algorithm with the linking length set to be 0.2 times the 
mean particle separation (Davis et al. 1985). The Hier- 
archical Bound- Tracing (HBT) algorithm recently devel- 
oped by Han et al. (2011) is then applied to split the 
FoF halos into disjoint, self-bound subhalos. The HBT 
algorithm finds and traces dark matter subhalos in sim- 
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ulations based on their merger hierarchy. Thus, when 
compared to previously existing subhalo finders, this new 
algorithm has the unique advantage that it is able to well 
resolve the subhalos in high density environment, while 
keeping strict physical track of their merger history. In- 
side each FoF halo, the subhalos are divided into two 
classes which are treated differently during tracing: 'cen- 
tral' subhalo and 'satellite' subhalos. The central subhalo 
is defined as the most massive subhalo of an FoF halo, is 
able to grow through accreting mass within the host halo 
and normally contains most of its mass. All the subhalos 
in an FoF halo except the central one are called satellite 
subhalos. 

We determine the position and velocity of each subhalo 
to be the center of mass and bulk velocity of its core, 
consisting of 25% of its bound particles with the lowest 
potential. For each central subhalo, we define the virial 
radius, i?200> as the radius within which the mean density 
of dark matter particles is equal to 200 times the critical 
value. The virial mass of the central subhalo is then 
defined as the mass within -R200: 

47T 

M 200 = 200 (Oc „ t — Rl 00 = lOOG-^^iifoo. (1) 

For a satellite subhalo, we use the A/200 a t the epoch 
when it was last the central subhalo in its host halo. 
Hereafter, we refer to both central and satellite subhalos 
as 'halos', and we refer to Af2oo defined in this way as 
the 'halo mass', if not specifically pointed out. 

3. RESULTS 

3.1. Group-galaxy cross-correlation functions 

In this work we estimate the velocity dispersion of 
galaxies by modelling the redshift distortions in group- 
galaxy cross-correlation functions. Thus, for each sub- 
sample of groups listed in Table 1, we begin by estimating 
the redshift-space cross-correlation function, ^ s '(r p ,7r), 
between the group central galaxies and the galaxies in the 
reference sample, using the commonly-used estimator 



&\r p ,ir) 



N R CG(r p ,w) 
N G CR{r p ,ir) 



- 1. 



(2) 



Here r p and ir are the separations perpendicular and par- 
allel to the line of sight; Nq and Nr are the number of 
galaxies in the reference sample and in the random sam- 
ple; CG(r pi ir) and CR(r p ,n) are the cross pair counts 
between the group subsamplc and the reference sample, 
and between the group subsamplc and the random sam- 
ple respectively, both with perpendicular separations in 
the bins log r p ± | A log r p and with radial separations in 
the bins 7r± \An. To reduce sampling noise, the random 
sample is constructed with a number of particles that is 
10 times the number of galaxies in the reference sample, 
i.e., N R = 10 x N G . 

An estimate of the projected cross-correlation function, 
w p (r p ), is then obtained by integrating £^(j" p ,7r) over 
the line-of-sight separations tt. We take 



w p (r p ) 



^(r p ,ir)d7r = J2& ) (r p ,7r i )ATr i 



where we choose 7r max = 40/i~ 1 Mpc as the outer limit 
for the integration depth (in order to limit noise from 
distant uncorrelated regions) so that the summation for 
computing w p (r p ) runs from 7Ti = — 39.5/i _1 Mpc to 
7T80 = 39.5/i~ 1 Mpc, given that we use bins of width 
AiTi = 1 Mpc. We note that we have tried with differ- 
ent outer limits for the integration depth and the result- 
ing w p (r p ) changes little. We have corrected the effect 
of fiber collisions in the SDSS data using the method 
described in Li et al. (2006c). 

In Fig. 1, we show the projected cross-correlation func- 
tions determined in this way for some of our group sub- 
samples. The errors on the cross-correlation function 
measurements are estimated using the bootstrap resam- 
pling technique (Barrow et al. 1984). We generate 100 
bootstrap samples from the group subsample and com- 
pute the cross-correlation functions for each sample using 
the weighting scheme (but not the approximate formula) 
given by Mo et al. (1992). The errors are then given by 
the la scatter of the measurements among these boot- 
strap samples. 

In all the cases the projected cross-correlation func- 
tion can well be separated into two parts: a steeper in- 
ner part at separations below around 1 /i Mpc and a 
flat outer part at larger separations. In the language 
of the halo model (e.g. Cooray & Sheth 2002), the inner 
part is dominated by the 'one-halo' term where the pair 
counts are mostly galaxy pairs in the same halo (central- 
satellite pairs in the case of cross-correlation functions 
as in this study), and the outer part is dominated by 
the 'two-halo' term where galaxy pairs are mostly in 
separate halos (central-central plus central-satellite pairs 
between different groups). The separation where the 
transition between the one-halo and two-halo terms oc- 
curs increases monotonically from ~200 ft. _1 kpc for the 
lowest-mass systems in our sample with M* ~ lO lo Af , 
up to ~1 h M.pc for the most massive systems with 
A/* ~ 1O 12 M0. The one-halo term shows a systematic 
change in shape with the mass of central galaxies, with 
steeper slopes at lower masses and flatter slopes at higher 
masses. This implies that the distribution of satellite 
galaxies in low-mass halos is more concentrated than that 
in high-mass halos. Moreover, the amplitude of the two- 
halo term increases with increasing mass, reflecting the 
tight correlation between the stellar mass of the central 
galaxies and the mass of their host dark matter halos. In 
contrast, the slope of the two-halo term seems to vary lit- 
tle, reflecting the known fact that the bias in the galaxy 
distribution is related in a very simple way to the bias in 
the distribution of dark haloes. The features seen in the 
group-galaxy cross-correlation functions are all well con- 
sistent with what we found in Yang et al. (2005b), but 
with smaller error bars thanks to the improvement of the 
data. 

In the plane-parallel approximation, the projected 
cross-correlation function w p (r p ) is expected to be di- 
rectly related to the real-space cross-correlation function 
£,cg(r) by an Abel transform, 



w p (r p ) 



£c a [Jrl+y*)dy = 2 



dr. 



I r 2 _ r 2 



(3) 



(4) 

Considering the one- and two-halo term features seen in 
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Fig. 1. — The projected cross-correlation function between the central galaxies of groups and the reference galaxies, estimated for some 
of the group subsamples listed in Table 1. The symbols with error bars are measured using the SDSS DR7 group catalogue. The solid lines 
are the best-fits of Eq. (4) to the data. The dashed and dotted lines are the small-scale component modelled by an NFW profile and the 
large-scale component that is a biased version of the projected linear two-point correlation function of dark matter in a WMAP7 cosmology. 
The stellar mass range of the group central galaxies and the best-fitting values of the model parameters are indicated in each panel (note 
that the parameter r s is in units of h~ 1 Mpc). The virial radius R200 given by the estimated halo masses listed in Table 1 is indicated by 
an vertical arrow in every panel. See the text for details. 
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the w p (r p ) measurements as discussed above, as well as 
previous theoretical studies of the cross-correlation func- 
tion between halo centers and dark matter particles (e.g. 
Hayashi & White 2008), we model the real-space cross- 
correlation function £, cg (r) as a linear combination of two 
components as follows, 

tcg(r) = £nfwM + b£ dm (r). (5) 

The first component, £nfw(0j is the NFW profile, 

Sr. 



£nfw(V) = 



(r/r s )(l + r/r s y 



where 



Sr.= 



200 



3fi m ln(l + c) - c/(l + c) 



(6) 
(7) 



The second component, £d m (V), is the linear two-point 
correlation function of dark matter at z — for the 
WMAP7 cosmology, obtained by Fourier transforming 
Pl(k), the linear power spectrum calculated using the 
CAMB code of Lewis et al. (2000). Thus, the model has 
three free parameters: c, r s and &, and we will call them 
concentration parameter, characteristic radius and bias 
factor in what follows, though they may have different 
meanings from their original definitions in the literature. 

For a given set of parameters, the model predicts 
both the real-space correlation function £ C n(?") through 
Eqn. (5) and the projected function w p (r p ) through 
Eqn. (4). We determine a best- fitting model for each of 
our group subsamples by comparing the predicted w p (r p ) 
with the observed one, and we define the best-fitting 
model to be the one giving a minimum \ 2 - O ur best 
models are plotted in Fig. 1 as solid black lines, and 
the two components of each model are plotted in red 
and blue dashed lines separately. The model parameters 
are also indicated. As can be seen, the two-component 
model provides a compact and accurate description of 
the observed w p (r p ) for all our group subsamples. It is 
interesting that all the model parameters show system- 
atic trends when one goes from the lowest-mass groups 
to the highest-mass groups, in the sense that both r s and 
b increase with mass, while the concentration parameter, 
c, decreases with mass. This is again consistent with the 
picture that more massive galaxies are hosted by more 
massive dark halos within which the galaxy distribution 
is less concentrated and extends to a larger radius. 

Moreover, the concentration parameter of satellite dis- 
tribution as inferred from our best-fitting models is 
smaller than the concentration of dark matter halos given 
by high-resolution simulations (e.g. Zhao et al. 2009), 
implying that the distribution of galaxies in dark halos is 
more extended than dark matter particles. This is consis- 
tent with some previous studies (e.g. Yang et al. 2005b; 
Chen et al. 2006; More et al. 2009a), though other stud- 
ies produced opposite results, especially for luminous 
red galaxies (e.g. Masjedi et al. 2006; Watson et al. 2010, 
2012; Tal et al. 2012). Nevertheless, we would like to em- 
phasize that our model parameters should not be over- 
interpreted, although they exhibit reasonable values and 
interesting trends with galaxy mass that seem to be eas- 
ily understood from our knowledge of galaxy-halo con- 
nections. 

In Fig. 2 we plot again the measurements and the best- 
fits of the projected cross-correlation function, as func- 



tions of both r p (left panel) and r p /i?200 (right panel). 
The systematic trends of w p (r p ) with central galaxy mass 
as seen from Fig. 1 are more clearly seen here. These in- 
clude the increase in the amplitude of the two-halo term 
on large scales, the change in the shape of the one-halo 
term on small scales, and the shift in the transition scale 
between the two terms on intermediate scales. 

3.2. Velocity dispersion profiles 

In the previous subsection we have obtained the group- 
galaxy cross-correlation functions, both in redshift space, 
(r p , 7r) and in real space, £ C g(?*)- In Fig. 3 we show the 
contours of £^ (r p , tt) for some of our group subsamples. 
Both the 'Finger-of-God' (FoG) effect on small scales and 
the infall effect (Kaiser 1987) on large scales are clearly 
visible: on small scales £,^(r p ,n) is stretched in the tt- 
direction and on large scales the contours are squashed. 
The FoG effect is caused by the peculiar, virialized mo- 
tions of galaxies within dark matter haloes, and is much 
more pronounced in the more massive haloes, reflecting 
their larger velocity dispersions. Here, we model these 
redshift-space distortions in detail, in order to infer the 
velocity field in and around dark haloes. 

Our method relies on the fact that the peculiar mo- 
tions of galaxies affect only their radial distances in red- 
shift space. Thus the information for peculiar velocities 
along the line of sight can be recovered by modelling the 
redshift-space cross-correlation function £,^(r p ,n) as a 
convolution of the real-space function £ C g(?") with the dis- 
tribution function of the peculiar velocity f(v cs ): 



Z {s) (r p ,Tr) = / f(v ca )t cg (yjri + in-vc)*) 



dVr 



(8) 



where v cs is the peculiar velocity of satellite galaxies rel- 
ative to the central galaxy. We adopt a Gaussian form 
for f(v cs ): 



f(Vcs) 



1 



exp 



■(Vcs ~ V C s) 



2 1 



(9) 



where is the mean and a v is the dispersion of the 
one dimensional peculiar velocities. Note that we have 
assumed that the distributions of the peculiar velocity 
and the velocity dispersion are isotropic. Assuming an 
infall model for t^j(r) , the velocity dispersion a v can then 
be estimated as a function of the projected separation r p 

by comparing the observed ^ s ' l ( r pj 7r ) with the modelled 
one. 

Following Croft et al. (1999), we use the non-linear 
spherical collapse model as our infall model, 

SJrY 

Sr.ut, 



v c ^(r) = v l {r) [1 + S(r)} 



-0.25 



exp 



(10) 



The first factor at the right-hand side, vi, is the infall 
velocity in the linear perturbation case, which depends 
on the distance to the local density maximum and is 
directly related to the density contrast (Peebles 1980): 



vi(r) 



1 



H n™rd(r), 



(11) 



where 6(r) is the overdensity inside radius r 
3J c 3 »(r) 



S(r) 



■TO 



£ cp (x)x z dx. (12) 
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Fig. 2. — Left: the Wp(r p ) measurements from the SDSS/DR7 group catalogue are plotted in symbols with error bars for some of our 
group subsamples, and are compared to the best-fits of Eqn.(4) to the data which are plotted in solid lines. Right: same as in the left panel 
except that the projected separation r v is scaled by the virial radius -R200 on the x-axis. The stellar mass ranges of the central galaxies are 
indicated. 



Here £ cp is the linear cross-correlation between galaxy 
clusters (or groups in our case) and matter. We use the 
simple linear biasing picture, so that 



jc P{r) = jm 



1 



£ cg (x)x 2 dx, 



(13) 



where £ C9 is the group-galaxy cross-correlation function 
in real-space and b is the linear bias factor in Eqn.(5), 
both of which have been obtained in the previous sub- 
section for each of our group subsamples. The power- 
law part at the right-hand side in Eqn.(lO) is intro- 
duced by Yahil (1985) (also used by Lilje & Efstathiou 
1989) as a good approximation to the exact solution 
for the non-linear collapse of spherically symmetric clus- 
ters. The exponential truncation in Eqn.(lO) is suggested 
by Croft et al. (1999) in order to take into account the 
highly non-linear behavior of the infall velocity at high 
overdensities. Tests against N-body simulations showed 
that Eqn.(lO) with S cu t = 50 provides an accurate de- 
scription of the average infall velocity of matter into 
galaxy clusters (Croft et al. 1999). 

In Fig. 3 we plot the contours of our best-fitting models 
of £^( r p> 7r ) obtained by applying Eqn.(8) to the data. 
In Fig. 4 we plot the velocity dispersion, a v , obtained by 
applying Eqn.(8) to the data, as functions of both the 
projected separation (r p , left panel) and the separation 
scaled by the virial radius (r p /i?200; right panel). There 
are several interesting trends/features that can be read 
from this figure. First, at fixed radii the velocity disper- 
sion is a strongly increasing function of central galaxy 



mass (thus halo mass), and the mass dependence be- 
comes stronger for more massive systems. The velocity 
dispersion within the scales of the one-halo term ranges 
from ~ 100 kms -1 for the lowest-mass systems in our 
sample (M* ~ 2 x 10 10 M o ), up to ~ 700 kms" 1 for 
the most massive systems with M* ~ 5 x 10 11 M Q . Sec- 
ond, the velocity dispersion shows a roughly flat profile 
at radii smaller than a few 100 /i _1 kpc, with a slight in- 
crease at the smallest radii for the massive systems with 
M* > 3 x lO n Af . This echoes the early results re- 
vealed from CNOC1 survey by Carlberg et al. (1997b) 
who found that for their clusters the line-of-sight veloc- 
ity showed a flat profile at radii larger than 0.3 virial 
radius with a very slight decline (see also Carlberg et al. 
1996; Biviano & Katgert 2004), as well as later studies 
on galaxy groups by Carlberg et al. (2001) who found 
flat or slowly rising profiles for small groups. Finally, 
the velocity dispersion profile exhibits a broad bump at 
around 1 Mpc, which is more pronounced for less massive 
systems. As we will show below, this is likely caused by 
contamination from neighboring larger halos where the 
high-speed galaxies contribute significantly to the veloc- 
ity dispersion measurement of the smaller system. 

3.3. Understanding the origin of the bump at around 1 

Mpc 

In order to understand the origin of the bump at 
around 1 Mpc in the velocity dispersion profile of the 
low-mass systems, we have performed the same analyses 
using our simulations and compared the results with the 
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Fig. 3. — The redshift-space group-galaxy cross-correlation function £( s ' (r p , tt) for groups with central galaxies in different stellar mass 
intervals, as indicated in each panel. The gray shaded regions show the measurements from the SDSS/DR7 data and the dashed red lines 
are the best-fits of Eqn.(8) to the data. 
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observation obtained above. 

We take all the halos from the z — outputs of 
our simulations, and we assign a stellar mass to each 
halo by applying the commonly- applied subhalo abun- 
dance matching model (SHAM; Vale & Ostriker 2004; 
Conroy et al. 2006; Shankar et al. 2006; Conroy et al. 
2007a; Baldry et al. 2008; Moster et al. 2010; Guo et al. 
2010; Neistein et al. 2011b,a; Yang et al. 2012). The 
model assumes that each halo has a galaxy at its center 
and the stellar mass of a galaxy is an increasing function 
of the maximum mass ever attained by its halo. In prac- 
tice, we obtain the relationship between dark halo mass 
-W200, as defined in § 2.3, and the stellar mass M* of the 
central galaxy by matching the number density of dark 
matter halos with mass above Af2oo with the number 
density of galaxies with stellar mass above M*. For this 
we have used the stellar mass function of galaxies in the 
local Universe measured by Li & White (2009) from the 
SDSS/DR7 galaxy sample, corrected in Guo et al. (2010) 
using "model" magnitudes. For the analysis in this sub- 
section we assume there is no scatter in the M200 — M* 
relationship, and we will take into account the effect of 
the scatter in later analyses. 

We note that, as pointed out recently by Yang et al. 
(2012), one problem of the current SHAMs is that they 
generally neglect the evolution of the central galaxy-host 
halo relation (see also Neistein et al. 2011a) and the dif- 
ferent stripping/disruption of satellite galaxies with re- 
spect to the subhalos. Realistically, satellite galaxies 
were accreted at higher redshifts and should have ex- 
perienced different evolutionary processes from the cen- 
trals. Nevertheless, as shown in Yang et al. (2012), cur- 
rent high redshift observations do imply that there is lit- 
tle evolution of the central-host halo relation at low red- 
shifts (z < 1). Moreover, the central- host halo relation 
for halos with masses > 10 12 M Q is found to be insensitive 
to these caveats. Thus the approach of SHAMs is quite 
safe for our study where the galaxies are at z ~ and 
are hosted by dark halos with mass well above 10 12 A/ Q 
(see below). 

The galaxy groups used here are identified using the 
so called halo-based group finder, where galaxies are 
grouped according to their common dark matter halos 
following the NFW profiles, i.e., by assuming the halos 
to be Spherical Ovcrdcnsity (SO) structures. Such kind 
of SO method can be quite different from the FoF algo- 
rithm where each halo/group is a self-bound virialized 
system. On the other hand, the FoF halos might have 
bridging effect — in some cases, two or more subhalos 
that are well separated in space are associated with a 
common FoF halo and only one of them (usually the 
most massive one) is classified as the central subhalo. In 
order to roughly take into account the difference between 
the FoF and SO halo finding algorithm, we have found 
out all the satellite subhalos that are located well be- 
yond the virial radius of their central subhalo (i?200,c), 
by requiring that d cs > i?200,c + R200,s, where d cs is the 
distance from the central subhalo and i?200,s is the virial 
radius of the satellite subhalo at the epoch when it was 
last the central object of its host halo (see §2.3). We re- 
identify such satellite subhalos as central subhalos and 
include them into the central galaxy subsamples when 
performing the analysis below. 



We have estimated both the cross-correlation func- 
tions and velocity dispersion profiles for the halo cat- 
alog. First, the central subhalos arc divided into a 
set of subsamples according to stellar mass, with the 
same mass intervals as listed in Table 1 for the real 
subsamples. Next, each subsample is cross-correlated 
with the full catalog to compute a redshift-space cross- 
correlation function, £( s )(r p ,7r), as well as a projected 
cross-correlation function, w p (r p ). We have followed the 
distant-observer approximation, defining the x — y plane 
of the simulation box as the 'sky' and the z-axis as the 
line of sight. The difference in the z-axis position be- 
tween two galaxies and the difference in their peculiar 
velocity along the same axis combine to give a line-of- 
sight separation, tt, thus including the redshift distor- 
tions in the computation of ^ s >(r p ,n). We should point 
out that we have added a random component to the pe- 
culiar velocity of the central subhalos, assuming it to 
follow a Gaussian distribution function with a la width 
of O.25V200, where V200 is the virial velocity of the halo. 
This takes into account the fact that the central galaxies 
of dark halos are not strictly at rest, but move relative 
to the halo center with a velocity that is 20-30% the 
virial velocity (Diafcrio et al. 1999; Berlind et al. 2003; 
Yoshikawa et al. 2003; van den Bosch et al. 2005; see 
also Skibba et al. 2011 and references therein). Finally, 
the velocity dispersion profiles are determined by apply- 
ing the ^ s >(r p ,ir) model in Eqn. (8) to the (r p , tt) 
measured from the halo catalogs, in exactly the same 
way as adopted for the real samples. 

In Fig. 5 we plot in blue dash-dotted lines the veloc- 
ity dispersion profiles measured from one of our simula- 
tions, the one with ag = 0.85, and compare these with 
measurements from the SDSS data (symbols with error 
bars). At all masses, the model and the data show very 
similar profiles over the full range of r p probed. In par- 
ticular, the bump at around 1 Mpc is seen clearly in the 
model. For systems with intermediate to high masses 
(M* > 5 x 1O 1O M0), the bump in the model is compa- 
rable to that observed, in both amplitude and shape. 
At lower masses, the bump is stronger in the model. 
As pointed out above, one possible reason behind the 
bump in the observed profiles is the contamination from 
neighboring halos in which the peculiar motions of galax- 
ies can contribute significantly to the velocity dispersion 
measurement of the system in question. To test this hy- 
pothesis, we have done a simple experiment in which we 
exclude all the halos from our central galaxy subsam- 
ples that are close to a neighboring larger halo, with the 
distance between the two halos smaller than the sum of 
their virial radii. The resulting velocity dispersion pro- 
files are plotted in red solid lines in Fig. 5. As can be 
seen, the bumps at around 1 Mpc become very weak or 
disappear, as expected. 

The change in the velocity dispersion profile is 
strongest at the lowest masses where the bump is no 
longer significant after the close halos are 'blended'. We 
have played around with different distance limits for 
blending the halos, and found that the bump becomes 
stronger for smaller blending radii. This can be under- 
stood from the fact that a smaller blending radius will 
include more satellite subhalos located in the inner re- 
gion of FoF halos where the contamination in velocity is 
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more serious. At the highest masses (M* > 2 x 10 11 M Q ), 
the velocity dispersion profile is essentially unchanged, 
reflecting the fact that the majority of the high-mass 
galaxies in the local Universe are central galaxies of dark 
matter halos. At these masses, the satellite fraction is 
found to be as low as ~ 10% by galaxy-galaxy lensing 
analyses of the SDSS data (Mandelbaum et al. 2006), as 
well as the SDSS-bascd halo occupation models (Cooray 
2006; Tinker 2007; van den Bosch et al. 2007). 

It is encouraging that the 'deblending' effect doesn't 
change the velocity dispersion profile within the virial ra- 
dius for all the stellar mass subsamples, except the lowest 
ones with log M*/Mq < 10.3. This shows that, although 
the profile beyond the virial radius is somewhat uncertain 
due to the contamination from neighboring larger sys- 
tems, the inner part can be well determined with the cur- 
rent data and methodology Therefore, in what follows 
we limit our analyses to masses above log M*/Mq = 10.3 
and focus on velocity dispersions within the virial radius, 
and we use the original subhalos identified by our HBT 
subhalo finder, without further mimicking the group find- 
ing algorithm of Yang et al. (2007). 

3.4. Constraining density fluctuation parameter as 
and the stellar mass vs. halo mass relation 

The velocity dispersion of galaxies reflects the action of 
the local gravitational field, and so should provide inde- 
pendent constraints on the connection between galaxies 
and dark matter halos. In addition, the relationship be- 
tween galaxy velocity dispersion and mass is expected to 
depend on the density fluctuation parameter as- This 
is simply because as largely determines the abundance 
of the most massive halos, thus the velocity dispersion of 
the most massive galaxies that are hosted by these halos. 
Here, we extend the analysis in the previous subsection 
and show that the as parameter and the stellar mass vs. 
halo mass relation can be simultaneously constrained us- 
ing our estimates of velocity dispersion obtained above. 

To the end we use all the simulations with different 
as values. We also allow the central galaxy mass of ha- 
los of given dark matter mass to have a certain amount 
of dispersion. Following the standard practice, we as- 
sume the dispersion in logM* to be Gaussian and in- 
dependent of halo mass, with rms values ranging from 
Cms = to 0.2 dex. The rms values exceeding about 
0.2 dex are excluded because they are inconsistent with 
the steep high-mass falloff of the stellar mass function of 
Li & White (2009). 

In Fig. 6 we plot the velocity dispersion profiles esti- 
mated from the subhalo catalogs for four of the stellar 
mass intervals (panels from left to right), and for mod- 
els with different dispersions in the stellar mass vs. halo 
mass relation (from top to bottom) and with different as 
values (the colorful lines). The observational results are 
plotted in symbols with error bars for comparison. The 
virial radii are indicated in each panel with the verti- 
cal lines, which are given by Eqn.(l) for the median halo 
mass of subhalos in a given stellar mass subsample. Over- 
all, the velocity dispersions within the virial radii depend 
on both as and a ms , but the behavior differs with stellar 
mass. At low masses, the velocity dispersion increases 
with increasing a ms , but shows very little dependence 
on as- At high masses, the velocity dispersion depends 
on both parameters with a stronger trend with as- The 



velocity dispersion is higher for systems with fixed mass 
and with fixed a ms if as is larger, well consistent with 
the naive expectation as discussed above. The weak de- 
pendence on as at the low-mass end is very encouraging, 
in the sense that the degeneracy between as and a ms as 
seen at high masses can be effectively broken, allowing 
the two parameters to be constrained simultaneously. 

For each stellar mass interval, we determined an aver- 
age velocity dispersion using the measurements falling in 
the range 0.3i?2oo < r p < ^?200- We exclude the inner- 
most part within 0.3i?2oo considering the facts that the 
velocity dispersion profile slightly increases at the small- 
est radii for massive systems (see Fig. 4) and this is not 
seen in the models at similar masses, and that we don't 
have velocity dispersion measurements on scales below 
0.3i?2oo for low mass systems (again see Fig. 4). The 
former fact is consistent with previous findings that the 
inner slope determined from gravitational lensing analy- 
ses is steeper than that expected from the NFW profile, 
implying that baryonic processes have modified the inner 
profiles (e.g. Koopmans et al. 2006). The average veloc- 
ity dispersions obtained from the SDSS data are listed 
in Column 5 of Table 1. In Fig. 7 we show the average 
velocity dispersion a v as a function of stellar mass. Re- 
sults from the SDSS data are plotted in symbols with 
error bars and are repeated in every panel. These are 
compared with results from the various models which 
have different values of as (panels from left to right) 
and a ms (different lines in each panel). The degener- 
acy between as and a ms , as well as the capability of the 
a v — M* relation for breaking the degeneracy, are seen 
more clearly in this figure. At low masses, the velocity 
dispersion shows strong dependence on a ms . The depen- 
dence on a ms is much weaker at high masses, particularly 
at ~ 2 x 10 n M Q where the velocity dispersion is essen- 
tially determined by as alone. The data at these masses 
clearly suggest that an intermediate value of as = 0.85 
works better than either higher or lower values. 

In order to have a more quantitative constraint on both 
as and a ms , we have performed the analysis for more 
models by varying the two parameters, with a ms ranging 
from to 0.2 dex with a step size of 0.01, and as ranging 
from 0.75 to 0.95 with a step size of 0.01. For given a ms , 
we estimate velocity dispersions for models of different 
as by interpolating the a v — M* relation obtained from 
the three simulations, assuming the velocity dispersion 
at fixed stellar mass to scale linearly with as- This as- 
sumption is reasonably true considering the following two 
facts. First, the change in velocity dispersion with as is 
fairly small, as can be seen from Fig. 7. Second, when ob- 
taining the number density of galaxy clusters through the 
Press & Schechter (1974) formalism, Mo et al. (2010, see 
their Eqns. 7.75 and 7.76) found that as scales with clus- 
ter mass as a 8 oc M^ 2 0^ /3 , where /3 w 0.6 + 0.8M7 m 
for CDM-type power spectra with index n = 1, leading 
to a nearly linear relation between as and velocity dis- 
persion a v for the assumed cosmology, a v oc erg 9 , given 
that the mass is roughly proportional to a% . 

We compare the velocity dispersions of each model 
to the SDSS measurements, and define the best-fitting 
model to be the one giving a minimum \ 2 computed as 
follows: 

X 2 (<7 8 ,cr ms )=X T C- 1 X. (14) 
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Fig. 6. — Velocity dispersion profiles estimated from subhalo catalogs for different stellar mass intervals (panels from left to right), and 
for models with different dispersions in the stellar mass vs. halo mass relation (panels from top to bottom) and with different density 
fluctuation parameter: erg = 0.75 (red solid line), 0.85 (green dashed line) and 0.95 (blue dash-dotted line). The estimates from the SDSS 
data are plotted in symbols with error bars for comparison. The virial radii from the models are indicated by the vertical lines. 
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Fig. 7. — Velocity dispersion averaged over 0.3i?200 < fp < ^200 is plotted as a function of central galaxy mass. Results from the SDSS 
data are plotted in symbols with error bars and are repeated in every panel. Results from the models are plotted in different lines, which 
have different values of erg (panels from left to right) and a ms (different lines in each panel). 
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Fig. 8. — x 2 > defined by Eqn.(14) and derived by comparing the 
average velocity dispersion within the virial radius of dark matter 
halos as predicted by the subhalo abundance matching model and 
as measured from groups of different central galaxy masses from 
the SDSS data, is plotted with respect to the minimum value Xmin 
(indicated by an 'x'), in the grid of the two model parameters, o~a 
and <T m s- The 68.3%, 90% and 95.4% confidence levels are plotted 
as red dashed lines, and the area between every two neighboring 
levels is filled with shaded regions. 
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Fig. 9. — Median halo mass M200 for central subhalos with given 
stellar mass is plotted as a function of the velocity dispersion mea- 
sured using our methodology, for the models with different values 
of erg and cr ms . The black solid line is the relation obtained by fit- 
ting to the result from the model with erg = 0.85 and <r ms = 0.15 
(plotted as big stars), and is adopted as our calibration for esti- 
mating halo masses from velocity dispersion measurements. The 
calibrated estimator is indicated at the bottom of the figure. 



Here, X = {Xj}(j = 1, m) is an to x 1 vector with 



A', 



■jdel 



SDSS 
J v,j 7 



(15) 



where to is the number of stellar mass subsamples, 
<f 55 is the velocity dispersion measured for the jth 
subsample from the SDSS and a 



■model 



is the result for 



the jth subsample from the model. The to x m matrix 



{Cij}{i,j 



1. 



, to) is the covariancc matrix of 



the measurements from the 100 bootstrap samples con- 
structed from the SDSS data, given by 



Cij — 



1 
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SDSS\ 
v,i ) 



boot 
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SDSS\ 



'•■J 



. (16) 



where n = 100 is the number of bootstrap samples, and 
(Ty°%l is the measurement for the ith stellar mass interval 
from the fcth bootstrap sample. 

Fig. 8 plots the contours of A% 2 = X 2 ~ Xmin m the 
grid of the two parameters. The 1, 2 and 3cr confidence 
regions, computed for to = 14 and 2 parameters, are in- 
dicated using the dashed lines plus shaded regions. The 
minimum Xmin appears at <rg = 0.86 and u ms = 0.16 
with x 2 Jd.o.f = 0.4, indicating that the fit is acceptable. 
There is a mild degeneracy between the two parame- 
ters in the sense that the models with smaller er 8 and 
larger er ms and the models with larger ag and smaller 
cr ms can both provide a reasonable fit to the data. While 
the minimum is the preferred solution at the la level, 
there is an elliptical degeneracy region extending to (0.82, 



0.17) and (0.90, 0.14), that are allowed at 2 - 3cr level. 
Our constraint on a ms is in good agreement with previ- 
ous studies which also found that the dispersion in stel- 
lar mass at fixed halo mass is ~ 0.13 — 0.20 dex (e.g. 
Wang et al. 2006, 2007; Yang et al. 2009; Behroozi et al. 
2010; Moster et al. 2010; More et al. 2011; Yang et al. 
2012). 

Our constraint on eg is slightly larger than the re- 
sult from the Wilkinson Microwave Anisotropy Probe 7 
(WMAP7) and low redshift supernova and baryon acous- 
tic oscillation data (Komatsu et al. 2011), which was 
(is = 0.816. We note that this difference can be largely 
explained by the slightly different matter density param- 
eter, which was Vl m = 0.275 in Komatsu et al. (2011) 
and is VL m = 0.268 in our simulations. This implies that 
our preferred value would be shifted to a% ps 0.84 for the 
WMAP7 matter density parameter, thus within the la 
level of the WMPAP7 result. 

3.5. Estimating galaxy halo masses 

The velocity dispersion is caused by the local gravita- 
tional field and so provides a direct measure of the dark 
halo mass of galaxies. We use our simulations to cal- 
ibrate the relationship between the velocity dispersion 
measured using our methodology and the halo mass for 
central subhalos of different stellar masses. As plotted in 
Fig. 9, the best-fitting relation is 

log(M 2O o/M ) = 7.O + 2.631og(cr l ,/kms" 1 ) 

-56.0/(cr t ,/kms~ 1 ), (17) 
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Fig. 10. — Dark matter halo mass as function of galaxy luminosity (left-hand panel) and stellar mass (right-hand panel). Black circles 
and red triangles are estimated from the velocity dispersion measurements, for all the central galaxies in our group catalog and for the 
subset of early-type centrals, respectively. Blue squares are the result obtained by Mandclbaum et al. (2006) by stacking the gravitational 
lensing signals of the early-type central galaxies in an earlier SDSS data release. 



obtained using the subhalo abundance matching model 
applied to the simulation with ct 8 = 0.85 and with the 
dispersion parameter a ms = 0.16. As shown above, this 
model provides the best match to the SDSS data in the 
stellar mass versus velocity dispersion relation. 

We apply the halo mass estimator in Eqn.(17) to 
our measurements of velocity dispersions. The result- 
ing halo masses are listed in the last column of Table 
1 and are plotted in Fig. 10 as a function of both lu- 
minosity (left-hand panel) and stellar mass (right-hand 
panel) . For comparison, we also plot the results obtained 
by Mandclbaum et al. (2006) from stacking the gravita- 
tional lensing signals from the SDSS data. We note that 
about 90% of our central galaxies are early-type galax- 
ies, which arc defined in practice to have the parameter 
f rac_deV > 0.5 from the SDSS photo pipeline. In or- 
der to make a fair comparison, we have measured the 
velocity dispersions for groups with early-type centrals 
and plot the results as red triangles in Fig. 10. The 
result doesn't change much when compared to that of 
the full sample. The halo masses of early-type cen- 
trals presented in Mandelbaum et al. (2006) are plotted 
as open squares. As can be seen, the halo masses esti- 
mated from galaxy velocity dispersions agree pretty well 
with those from gravitational lensing signals, over the 
full luminosity range (2 x 1O 1O L < L < 1.5 x 10 n M o ) 
where the halo mass estimates are available from both 
studies. Note that we have shifted the halo masses of 
Mandclbaum et al. (2006) by about -0.1 dex in order 
to take into account the different halo mass definitions 
adopted in the two studies. 

The agreement becomes slightly worse when the 
halo mass is plotted 8jS 0, function of stellar mass, 
as shown in the right-hand panel of Fig. 10. In 



this case our halo mass is systematically higher than 
that from Mandelbaum et al. (2006) by about 50% 
at intermediate to low masses. The stellar masses 
used in Mandclbaum et al. (2006) were estimated by 
Kauffmann et al. (2003) based on the SDSS spec- 
troscopy, while those in this study are based on the pho- 
tometric properties. Therefore, for this comparison in 
Fig. 10, we have shifted the former by about 0.1 dex 
in stellar mass; the correction factor was obtained in 
Appendix A of Li & White (2009) using a large sam- 
ple of galaxies for which both stellar mass estimates are 
available (see also fig. 17 of Blanton & Rowcis 2007). As 
can be seen in Li & White (2009), the difference be- 
tween the two stellar masses depends on stellar mass in 
a complicated way with varying scatter. Thus, a con- 
stant correction cannot fully get rid of the difference, and 
the larger discrepancy between our results and those of 
Mandelbaum et al. (2006) as expressed in terms of stellar 
mass must be partially (if not totally) due to this effect. 

4. SUMMARY AND DISCUSSION 

Using a large catalog of galaxy groups selected from the 
Sloan Digital Sky Survey (SDSS) data release 7 (DR7) by 
Yang et al. (2007), we have derived the velocity disper- 
sion profiles for groups with different masses of central 
galaxies, by modelling the rcdshift distortions in group- 
galaxy cross-correlation function. Our main results can 
be summarized as follows. 

• An Navarro et al. (1997) profile on small scales plus 
a biased version of the linear mass autocorrelation 
function can well describe the observed w p (r p ) over 
all the scales probed (15kpc < r p < 30Mpc), thus 
providing an accurate determination of the real- 
space cross-correlation function £ C g(f)- This has 
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been well expected from theoretical studies of halo- 
mass and galaxy-mass cross-correlation functions 
(e.g. Hayashi & White 2008). 

• The velocity dispersion within virial radius (-R200) 
shows a roughly flat profile, with a slight increase 
at the smallest radii (< 0.3i?2oo)- This is consis- 
tent with the theoretical expectation that baryonic 
condensation has modified the inner profile. 

• The average velocity dispersion within the virial 
radius, a v , is a strongly increasing function of cen- 
tral galaxy mass, ranging from ^130 kms -1 for 
the lowest-mass groups in our sample which have 
central galaxy masses of M* ~ 3 x 10 10 M Q , up 
to ~650 kms -1 for the highest-mass systems with 
M* ~ 5 X lO n M . The mass dependence is more 
remarkable at higher masses. 

We have extended our analysis of the SDSS data to 
a set of iV-body simulations with the concordance A 
cold dark matter (ACDM) cosmology but different val- 
ues of the density fluctuation parameter eg- We have 
assigned a stellar mass to each of our subhalos under 
the plausible hypothesis that the stellar mass of cen- 
tral galaxies is an increasing function of the dark matter 
mass of their host halos, usually called subhalo abun- 
dance matching model in the literature, and measured 
the velocity dispersion profile for the subhalos in dif- 
ferent stellar mass intervals. By comparing the results 
to measurements from the SDSS data, we have shown 
that the a v — M* relation obtained from the data pro- 
vides stringent constraints on both the a& parameter 
and a ms , the dispersion in log A/* of central galaxies 
at fixed halo mass. Our best-fitting model suggests 
that as = 0.86 ± 0.03 and a ms = 0.16 ± 0.03. The 
latter is in very good agreement with previous studies 
where a rns was inferred from semi-analytic catalogs of 
galaxy formation models (e.g. Wang ct al. 2006, 2007), 
from group catalogs (Yang et al. 2009), from satellite 
kinematics analysis (More et al. 2011), from clustering 
analysis (Mostcr et al. 2010), and from subhalo abun- 
dance matching (Behroozi et al. 2010; Guo et al. 2010; 
Yang et al. 2012). 

Our velocity dispersion measurements also provide a 
direct measure of the dark matter mass of the host halos 
for central galaxies of different luminosities and stellar 
masses, which are in good agreement with the results 
obtained by Mandelbaum et al. (2006) from stacking the 
gravitational lensing signals of the SDSS galaxies. The 
results also broadly agree with recent estimates by many 
other authors (e.g. van den Bosch ct al. 2004; Yang et al. 
2007; More et al. 2009b; Guo ct al. 2010; More et al. 
2011), though with varying degrees of subtle discrepan- 
cies due to different methodologies and date sets. 

It is encouraging to see that the velocity dispersion 
as a function of central galaxy mass as measured in our 
study is able to constrain the density fluctuation param- 
eter as, although the preferred value is slightly higher 
than that from the most recent result from WMAP7 
(Komatsu et al. 2011). We suspect that this difference 
can be explained if the effect of the different matter den- 
sity parameters adopted in the two studies is considered. 
It has been recently emphasized that the redshift-spacc 
distortion on very large scales (near or beyond the baryon 



acoustic oscillation scale), expected to be well detected 
in next-generation redshift survey like the BigBOSS ex- 
periment (Schlegcl ct al. 2011), will provide powerful 
constraints on the growth factor of structure and the 
physical properties of dark energy (e.g. Okumura & Jing 
2011; Samushia et al. 2012). Our study suggests that 
the redshift distortions on galaxy and group scales pro- 
vide complementary information useful for constraining 
cosmological parameters. When combined with more ob- 
servational measurements like the autocorrelation func- 
tion and pairwise velocity dispersion of galaxies as func- 
tions of their physical properties (e.g. Zehavi et al. 2005; 
Li et al. 2006a,b), our results should be able to provide 
more accurate and reliable constraints on as- We will 
come back to this point in future studies. 

There is still room to improve our study. First, we have 
assumed isotropic distributions for the peculiar velocity 
and velocity dispersion of satellite galaxies and this is 
roughly, but not exactly correct. Models incorporated 
with radial anisotropy might provide better determina- 
tions of the velocity dispersions. Using red, early-type 
galaxies instead of the whole population as the tracer 
may also help, as they are more likely consistent with 
isotropic orbits than are blue and late-type galaxies (e.g. 
Carlberg et al. 1997c; Biviano & Katgert 2004). Second, 
we have used groups with three or more member galaxies 
and assumed that our determinations of velocity disper- 
sion are insensitive to the richness limit of the group sam- 
ple. Larger and deeper surveys in future would allow us 
to test this assumption and improve our method. Third, 
previous studies found that Brightest Cluster Galaxies 
(BCG) may be not at the true center of their host halo 
(e.g. Skibba et al. 2011), and this effect could also intro- 
duce some sort of systematics. Fortunately, this is un- 
likely a significant issue for us, as such offsets seem to be 
small in general (a recent study by von der Linden et al. 
2012 found a median value of 20 kpc in the ffsets between 
X-ray centroids and BCGs in clusters). Fourth, we have 
assumed the dispersion in stellar mass is independent of 
halo mass and this might be tested and improved via 
semi-analytic galaxy formation models or halo occupa- 
tion distribution models. Finally, our work is limited to 
relatively large systems with M* > 3 x 10 10 M Q and could 
be extended to lower masses using deeper surveys. Pre- 
vious studies of the pairwise velocity dispersion (PVD) 
of galaxies revealed that the relative velocities between 
galaxies exhibit a well-defined minimum at intermediate 
luminosities/masses, and both fainter and brighter galax- 
ies have very high velocities, indicative of massive halos 
of cluster size (Jing & Borner 2004; Li et al. 2006a). It 
would be interesting to directly measure the halo mass 
for such faint galaxies by applying our methodology to 
groups of lower masses, e.g. down to - 6 x 1O 9 M 
using the Galaxy And Mass Assembly survey (GAMA; 
Driver et al. 2009). 
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